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Abstract 

A precise understanding of the radio emission from extensive air show- 
ers is of fundamental importance for the design of cosmic ray radio de- 
tectors as well as the analysis and interpretation of their data. In recent 
years, tremendous progress has been made in the understanding of the 
emission physics both in macroscopic and microscopic frameworks. A 
consistent picture has emerged: the emission stems mainly from time- 
varying transverse currents and a time- varying charge excess; in addition, 
Cherenkov-like compression of the emission due to the refractive index 
gradient in the atmosphere can lead to time-compression of the emitted 
pulses and thus high-frequency contributions in the signal. In this article, 
I discuss the evolution of the modelling in recent years, present the emis- 
sion physics as it is understood today, and conclude with a description 
and comparison of the models currently being actively developed. 



1 Introduction 

Radio emissi on from cosmic ray air showers has been measured for the first time 
in the 1960s ( Jellev et ah . 1965f) . This started a flurry of activities both on the 



experimental and theoretical side; for a review of this early history, please see 
(Allan, 1971). However, difficulties in the interpretation of the measurements 
and limitations of the analogue detection technique led to a complete cease of 
activities in the 1970s. 

The field was re vived in the early 200 0s by the LOPES (|Falcke et al.l . |2005|) 



and CODALEMA ()Ardouin et all . |2005|) experiments. It was clear from their 



start that a sound understanding of the radio emission physics was imperative 
for the interpretation of measurements and the optimization of the experimental 
designs. Therefore, in parallel to the revival of experimental activities, also 
theorists started new efforts for the modelling of air shower radio emission. 

Here, I give an overview of these "modern" efforts. I shortly review how 
the models have evolved and describe the approaches currently employed. Af- 
terwards, I give a summary of the current understanding of the radio emission 
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physics and conclude by presenting some results gathered by a task force aimed 
at comparing the existing models in detail. 



2 The evolution of modern radio emission mod- 
elling 

A driver for the modern efforts in radio de tection of cosmic ray a ir showers has 
been the paper by Falcke & Gorham [Falcke and Gorhaml . I20Q3I ). It predicted 
that air shower radio emission should be measurable with digital radio detectors, 
and presented a first estimate of the field strengths to be expected on the basis 
of the "geosynchrotron" approach. 



2.1 Obsolete approaches 

As a direct follow-up to this paper, the geosynchrotron co ncept was followed 



up w ith semi-analytic calculations in the frequency domain (IHuege and Falcke . 



2003 ) and with Monte Carlo calculations in the time domain (|Suprun et al 
20031 ) . In parallel, full-fledged Monte Carlo si mulation codes were developed, in - 



parti cular REASl (|Huege and Falckel l2005allb[ ) and ReAires (IDuVernois et al 



2005[) . REASl made the transition to REAS2 (|Huege et all . l2007l) by replac- 



ing the parameterized air shower model with one based on deta iled CQRSIKA- 
derived histograms. Other s emi-analytical calculations followed (jMever-Vernet et al 



2008t IChauvin et al.l . I2OI0I) . For a review of the evolution of these models in 



the time from 2003 to 2009, please see (iHuegel . [2OO9I) . 

However, it turned out that all of these models were employing an incomplete 
modelling of the radio emission physics. They did not take the radiation asso- 
ciated with the varying number of charges in air showers into account correctly, 
and therefore gave inconsistent results. A particular signature of this incon- 
sistency was that the models always predicted unipolar time-domain pulses, 
whereas other approaches predicted bipolar pulses. For a detailed discussion of 
the cause for and so lution of these discrepancies, I kindly refer the reader to 



(jHuege et al.l . l2012al ). Readers should be aware that the models listed in this 



section should be considered obsolete and should thus no longer be actively used. 
Still, many of the predictions that w ere made by these mo dels, e.g. on coher- 



ence e ffects (jHuege and Falckel 120031 ). air shower geometry (IHuege and Falcke . 
2005bl) . and energy and depth of shower maximum sensitivity (jHuege et al 
2008h were qualitatively correct. 



2.2 Current approaches 

With the understanding of the inconsistencies in the models, the field has made 
tremendous progress. Nowadays, many different approaches exist, and they 
all agree in their qualitative description of the radio emission features. The 
approaches can be grouped in the two categories of microscopic and macro- 
scopic models. The microscopic models follow individual particles (electrons 
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and positrons) in an air shower and calculate the associated radio emission. 
Macroscopic models rather describe the radio emission physics on the basis of 
macroscopic quantities such as currents and net charge. Both have their virtues 
and limitations. 



2.2.1 Microscopic approaches 



A number of microscopic models exist today. The formalisms used to calculate 
the electromagnetic radiation differ, but they have in common that no assump- 
tions are made on the actual "radio emission mechanism" . The electromagnetic 
radiation is calculated on the basis of classical electrodynamics applied in a 
specific formalism. The radiation that is emitted is thus directly governed by 
the movement, and in particular acceleration, of the individual particles. (In 
this sense, it might be more appropriate to call these approaches "simulations" 
rather than "models".) The following is a list of existing and actively main- 
tained microscopic models with a short description for each of them. (Other 
microscopic models exist but are no longer actively maintained and thus not 
listed here.) 



REA S3.1 (iLudwig and Huegel . l201ll) uses the endpoint formalism (jJames et al. 



20111 ) to calculate the radio emission on the basis of histogrammed COR- 



SIKA showers (the same as in REAS2) in the time-domain. 

SELFAS2 (jMarin and Revenul |2012|) calculates the radio emission on the 



basis of a formalism separating the electric field into components for the 
static contribution, charge variation and current variation in the time- 
domain. Its air shower model is based on the universality of histogrammed 
CORSIKA simulations (jlafebre et al.U2009l) . so no air shower simulation 



is needed for a radio simulation. 



• CoREAS (IHuege et al.ll2012bl ) employs the endpoint formalism as in REAS3.1, 
but calculates the radio emission directly in CORSIKA on a per-particle 
level, without any intermediate histogramming step. It also works in the 
time-domain. 



ZHA ireS ( Alvarez-Mufiiz et al. , 2012f) employs the ZHS algorithm ( Zas et al 



I1992 !) on a per-particle level in AIRES in either the time- or the frequency- 
domain. 

All of these models include a realistic treatment of the atmospheric refractive 
index. The two "full Monte Carlo" simulations CoREAS and ZHAireS describe 
the underlying air shower with the highest degree of complexity. 



2.2.2 Macroscopic approaches 

In contrast to the microscopic models, which calculate the emission from indi- 
vidual air shower particles and in principle make no assumptions on the emission 
"mechanism" , the macroscopic models calculate the emission as emanating from 
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the "bulk features" of the particle distributions, such as currents, net charge or 
dipole moments. Their advantage is that they are faster and give useful insights 
in the emission physics. On the other hand, in some cases free parameters have 
to be set, and it is sometimes difficult to separate the different contributions 
(| James et al.U201ll) . 



MGMR (jScholten et al.l . 120081 ) calculates the radio emission on the basis 
of time- varying transverse currents, time varying net charge and a time- 
varying dipole moment in the time-domain. The air shower model is based 
on parameterizations, some important parameters such as the drift veloc- 
ity of the particles and the shower-disk thickness need to be chosen. 

EVA (|Werner et al.U2012l ) employs the same time-domain emission calcu- 
lation as MGMR, but parameterizes the air shower particle distributions 
on the basis of individual CONEX simulations. It thus has no free param- 
eters to set. The parameterizations are currently one-dimensional, but 
could be generalized to multiple dimensions. 



• The semi-analytical model developed by Dave Seckel (|Seckei 120121 ) em- 
ploys a similar approach as MGMR, but works in the frequency domain. 



Of these models, 
refractive index. 



EVA also includes a realistic treatment of the atmospheric 



3 The emission physics 

Having shortly reviewed the models available and actively maintained today, it 
is instructive to discuss our current understanding of the radio emission physics. 

3.1 Identified contributions 

While disentangling different "mechanisms" for the radio emission from exten- 
sive air showers (or any complex ra diating system) is difficult and sometimes 



even misleading (jJames et al.l . 120111) , it is still often instructive to try. For air 



shower radio emission, the following contributions could be identified: 

• The dominant contribution to the radio emission is of geomagnetic origin. 
The emission is linearly polarized with the electric field vector aligned 
in the direction of the Lorentz force {v x B), independent of the ob- 
server location as shown in Figure [1^). The radiation can be under- 
stood as emanating from a time- varying transverse current. The current 
is a consequence of the equilibrium between acceleration of electrons and 
positrons in the geomagnetic field and deceleration by interactions with 
atmospheric molecules. The time-variation arises from the growth and 
decline of the number of charges during the air shower evolution. This is 
the mechanism that was already proposed in the 1960s by Kahn & Lerche 
( Kahn and Lerche . 19661 ). 



4 



• A second contribution to the radio signal arises from the time-variation 
of the net charge present in an air shower. This is the same mechanism 
responsible for radio emission of showers in dense media. In air, however, 
the contribution is only sub-dominant with respect to the geomagnetic 
emission. The emission is linearly polarized with the electric field vectors 
oriented radially with respect to the shower axis. Consequently, the polar- 
ization angle varies with observer location as presented in Figure [TJd). The 
net charge in an air shower arises mainly because electrons are knocked 
out of atmospheric molecules and are swept along with the shower disk. 
It is time-varying also because of the growth and decline of the number 
of charges in the air s howe r . Thi s is the emission mechanism proposed by 
Askaryan ( Askarvan . 1962 . 1965[ ). 



• To which extent the originally proposed "geosynchrotron" emission, i.e. 
the radiation associated directly with the acceleration of electrons and par- 
ticles in the geomagnetic field, contributes to the radio signal is not yet 
fully resolved; it is conceptually difficult to disentangle completely from 
the geomagnetic effects leading to the transverse currents. It might play a 
role at very high frequencies, where recent simulations have predicted po- 
larization characteristics s imilar to what has b een predicted in the earlier 
"geosynchrotron" models (jHuege et al.L l2012b[ l. 



• It has been known for a long time that atmospheric electric fields can in- 
fluence the radio emission fro m air shower s, and the effects have also been 
studied theoretically (Buiti nk et al. . 2010[ ). In recent modelling efforts, it 



could be shown that, depending on the relative orientation of electric field 
and shower axis, the transverse current and charge excess contributions 
can be amplified or dampened by atmospheric electric fields, see Figure [2j 

• An effect first noticed by Dave Seckel in REAS3 simulations is a secondary 
pulse arising at later times for suitable air shower geometries. This arises 
from the ground impact of the particles, which are stopped instantaneously 
when reaching the observation level in the REAS3 simulation. A more 
realistic treatment should take into account the electric properties of the 
soil and model the "transition radiation" arising from this ground impact. 
This could well be relevant at kHZ frequencies. 

• Some macroscopic models include other contributions such as a time- 
varying dipole or the ions which are left behind in the wake of the air 
shower. This might be seen as an illustration of the complexity of "sum- 
ming up mechanisms" . 

In summary, the main mechanisms for the radio emission are the time-varying 
transverse currents and the time-varying net charge present in the air shower. 
The superposition of these differently polarized contributions directly leads to 
prominent asymmetries in the radio emission footprint (in particular an east- 
west asymmetry for vertical air showers), as is obvious from Figure [1] 
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Figure 1: Alignment of the electric field vectors as a function of observer po- 
sition (shower core at the center) for a vertical air shower. Left: transverse 
current emission, right: charge-excess emission. Both contributions are linearly 
polarized. 
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Figure 2: North-south component of the electric field for an observer 250 m 
north of the core of a 30° inclined 10^® eV air shower simulated at the LOPES 
site with CoREAS 1.0 in the presence of various vertical atmospheric electric 
fields. 

3.2 Refractive index effects 

Another effect is important for the radio emission from extensive air show- 
ers. The refractive index of the atmosphere is not unity but has a value of 
~ 1.000292 at sea level and decreases with height as the atmospheric density 
decreases. The fact that the index is larger than unity means that for suitable 
geometries radiation emitted at different times and locations can reach a given 
observer at the same time. In such situations, the pulses are thus compressed in 
time, and consequently the frequency spectra of the emitted radiation can ex- 
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tend up to very hiRh frequencies well in the GHz regime (jAlvarez-Mufiiz et al 



2OI2I: iDe Vries et all . I2OIII: IWerner et all. 120121). Many such hiRh-frequency 
results have recently been presented (|Huege et all . l2012bl: [Carvalho Jr.l [2oTa 



De Vries et al.l . l2012f) . and they all agree qualitatively in their prediction of 
a "Cherenkov ring" appearing at high frequencies. The diameter of this ring 
is related to the geometrical distance between source and observer, and thus 
carries information on the air shower evolution. While it is adequate to label 
these refractive index effects "Cherenkov effects" , this should not be mixed up 



with what is usually referred to as classical "Cherenkov radiation" (jJames et al 
20111) . 



4 Where we stand today 

The modelling of radio emission from extensive air showers has come a long 
way, and with the resolution of the discrepancies present in the first "modern" 
approaches, a consistent picture has emerged. The goal of the modelers should 
be to understand the radio emission on a level on par with the systematic 
uncertainties of existing and future experiments, which means on a « 10% 
level. 

4.1 Model- model comparisons 

One way of judging how close the models are to reaching this goal is to compare 
them systematically. A task force involving authors of all models discussed in 
section [22] has started to undertake this effort, initiated at a workshop at Ohio 
State University in early 2012. In preparation for the ARENA conference, these 
comparisons were updated once more, and I have the honor of presenting a small 
excerpt of these on behalf of the task force in Figures [3l |4] and M 

In these figures, it is striking that the results from CoREAS and ZHAireS 
are very similar. These two codes model the interaction and emission physics 
with the highest complexity, i.e. on a microscopic per-particle level, and are in- 
dependent in their technical implementations of both the air shower simulation 
(CORSIKA vs. AIRES) and the radio emission calculation (endpoint formalism 
vs. ZHS algorithm). Their correspondence is thus a very promising indication 
that the models are indeed converging. REAS and SELFAS show slight de- 
viations from the full Monte Carlos, which might be related to the fact that 
they are based on information from histogrammed particle distributions rather 
than individual particles. Another striking observation is that the macroscopic 
MGMR and EVA models predict significantly higher amplitudes, especially for 
observers close to the shower axis. These differences to the microscopic models 
are not yet understood. 
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Figure 3: West-component of the electric field pulses for a vertical 10^^ eV 
proton-induced air shower at the site of the Pierre Auger Observatory as sim- 
ulated with the indicated models. The refractive index of the atmosphere was 
set to unity. 



4.2 Model-data comparisons 

While comparisons between different models are important and studying the 
reasons for deviations will be imperative to further improving the models, it is 
comparisons with data which matter in the end. Unfortunately, currently no 
data set is publicly available for the modelers to compare their models against. 
It would be of great use if experiments could publish a data set with all necessary 
details (including not only air shower parameters such as energy and angles but 
also the exact core positions and relative antenna locations) so that modelers 
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Figure 4: Absolute amplitudes of the frequency spectra for a vertical 10^^ eV 
proton-induced air shower at the site of the Pierre Auger Observatory as sim- 
ulated with the indicated models. The refractive index of the atmosphere was 
set to unity. 



can use them to benchmark their predic tions. Until then, co mparisons carried 



out by the experimental collaborations (jLudwig et al 
way to judge the models. 



20121 1 will be the only 



4.3 Open questions 

One issue that bothered modelers over the past year was the question whether 
the established calculations worked correctly very near the Cherenkov angle of 
individual particle tracks. A number of contributions to these proceedings have 
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Figure 5: West-component of the electric field pulses for a vertical lO^'' eV 
proton-induced air shower at the site of the Pierre Auger Observatory as sim- 
ulated with the indicated models. The refractive index of the atmosphere was 
modelled according to atmospheric density. 



addressed this issue and the preliminary conclusion is th at for th e geome tries 
relevant for air showers, the treatments seem to be valid ( James . 2012; B elovL 
Other aspects that should be studied in more detail than they have 



before include the transition radiation arising from the particle ground impact, 
the influence of realistically modelled atmospheric electric fields and possible 
reflection or scattering of radio waves on the particle plasma in the shower disk 
itself. 



5 Conclusions 

The modelling of radio emission from cosmic ray air showers has made tremen- 
dous progress in the past few years. A consistent picture has emerged, and 
several microscopic and macroscopic models are available and agree in their 
qualitative predictions. In the next few years, the models have to be bench- 
marked against each other and against measured data, with the goal to reach 
and prove a modelling accuracy within the systematic uncertainties of experi- 
ments. 
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